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X-RAY FADING AND EXPANSION IN THE “MINIATURE SUPERNOVA REMNANT” OF GK PERSEI 
D. Takei 1 , J. J. Drake 2 , H. Yamaguchi 3,4 , P. Slane 2 , Y. Uchiyama 5 , & S. Katsuda 6 

ABSTRACT 

We report on a second epoch of Chandra X-ray imaging spectroscopy of the spatially-resolved old 
nova remnant GK Persei. An ACIS-S3 observation of 97.4 ks was conducted in November 2013 after 
a lapse of 13.8 years from the last visit in 2000. The X-ray emitting nebula appeared more faint and 
patchy compared with the first epoch. The flux decline was particularly evident in fainter regions and 
the mean decline was 30-40% in the 0.5-1.2 keV energy band. A typical expansion of the brightest 
part of the remnant was l'/9, which corresponds to an expansion rate of 0'/14 yr^ 1 . The soft X-ray 
spectra extracted from both the 2000 and 2013 data can be explained by a non-equilibrium ionization 
collisional plasma model convolved with interstellar absorption, though do not allow us to constrain 
the origin of the flux evolution. The plasma temperature has not significantly evolved since the 2000 
epoch and we conclude that the fading of the X-ray emission is due largely to expansion. This implies 
that recent expansion has been into a lower density medium, a scenario that is qualitatively consistent 
with the structure of the circumstellar environment photographed soon after the initial explosion more 
than a century ago. Fainter areas are fading more quickly than brighter areas, indicating that they 
are fainter because of a lower ambient medium density and consequently more rapid expansion. 

Subject headings: stars: novae, cataclysmic variables — stars: individual (GK Persei, Nova Persei 
1901) — X-rays: stars 


1. INTRODUCTION 

Structures resulting from stellar explosions represent 
space laboratories that elucidate the physics of enor¬ 
mous blasts under extreme conditions, and are of fun¬ 
damental importance to astrophysics. The representa¬ 
tive of this kind of structure is a supernova remnant 
(SNR). X-ray observations of SNRs have provided wide- 
ranging insights into their plasma physics, s hock p he¬ 
nome n a, and into Cosmic-ray origins le.g.. iSla ne et all 
19991: lUchivama et al.l I2007t lYamaguchi et al.l 120091 : 
Katsuda et al.l 120131 ). However, the time history of an 
individual SNR is generally difficult to track because 
of the long timescales — comparable to the entirety 
of human written history — over which typical rem¬ 
nants evolve. With the notable exception o f specta cular- 
nearby r emn ants, s uch as Cas A (e.g., iFesen et~ahl 12 0011 
iPatnaude & Fesenl 12009( 1 , our understanding of remnant 
expansion and evolution has therefore largely depended 
on theory, numerical simulations, and systematic snap¬ 
shot observational assessments of objects sampling dif¬ 
ferent evolutionary phases. 

This paper presents a different approach to studying 
the X-ray evolution of an explosion remnant. We study 
not a supernova but a “classical” nova resulting from a 
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cataclysmic explosion in an accreting binary system com¬ 
prising a white dwarf and a red dwarf or giant compan¬ 
ion. The released energy and material propagate through 
the circumstellar environment similar to, but on much 
smaller scales than, those in supernovae. Classical no¬ 
vae are about 10 5 times smaller in energy, and typically 
factors of 10 5 -10 7 times smaller in ejecta mass. The anal¬ 
ogy with supernovae renders a nova remnant a miniature 
SNR with a short life: its full dynamical evolution occurs 
on a timescale comparable with a human lifetime and can 
potentially be followed in X-rays by multi-epoch imaging 
spectroscopy. 

Our target is Nova Persei 1901 (hereafter GK Per), 
which hosts the largest known X-ray emitting nebul a 
centered on a white dwarf binary dSeaauist et al.l H989al . 
In the following, we start by comparing Chandra Cycle-1 
data of GK Per obtained in 2000 with our recent observa¬ 
tion in Cycle-15 obtained in 2013. The earlier data have 
been analyzed in detail bv lBalmanl (12005D who character¬ 
ized the shock conditions, energetics of the remnant, and 
the shocked mass. Here, we use the second epoch data 
to examine the remnant expansion, cooling and fading, 
and discuss the results in the context of the physics of 
expanding explosion remnants. 

2. TARGET (GK PERSEI) 

GK Per is a mag netic w hite dwarf binary at the dis¬ 
tance of 470 pc dMcLaughlinl [i960) that un derwent a 
classi cal nova outburst on 1901 February 21 dWilliamsI 
HMD. Light echoes from the event photographed in 
1901 and 1902 illuminated a com plex of shell-like struc¬ 
tures arou nd the cent r al obj ect (lRltchevlll902t iPerrind 
Il902t) that iBode et all (119871 ) and iBode et al.1 (|2004H in- 
terpreted as part of an old planetary nebula. This neb¬ 
ula is also apparent as a vaguely bipolar-shaped struc¬ 
ture, elongated in the northwest-southeast direction, in 
the light of Ha, [O III] A 5007, and the IRAS 100 /im 
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band. iSchaeferl (1198811 noted that the lack of the sort 
of reflection features seen in GK Per around other novae 
indicates that the grain density around GK Per is unusu¬ 
ally high in keep ing with the planetary nebular origin 
suggested bv iBode et al.1 (119871 ). 

The extended nova remnant and its evolution has been 


imaged extensively 
(e.g., Seaouist et al. 

in the optical over several decades 
1989bt lAnupama & Prabhu 1993i 

Slavin et al] 1995); 

Lawrence et all 19951 Shara et al.1 

2012; Liimets et al.l 

2012j). It has a pronounced asyrn- 


metric appearance, with brighter emission to the south¬ 
west and a dearth of emission in the east. The limb 
defined by the extent of optical knots exhibits a decid¬ 
edly flattened appearance to the northwest and north¬ 
east. Hubble Space Telescope ( HST ) observations in 
1995 and 1997 revealed the details of the remn ants com- 
plex structure comprising a myriad of knots (IShara et al.l 
120121 1. Its present-day angular size based on optical im¬ 
ages from 2011 December 13 is about V in radius, and the 
kno ts are observed to be e xpanding with circular symme - 
try (iLiimets et al.ll2012h . lAnupama fe Kantharial (|2005ll 
found a possible deceleration of the ejecta after 1950 for 
the period 1901-2003, suggesting that the remnant is in 
a transition to an adiabatic phase. A typical expansion 
rate of the outermost optical knots in 2004-2011 was 0'/3— 
0 "5 yr - 1 , corresp onding to a velocity of 600-1000 km s _1 
(|Liimets et al.ll2012f). which i s simi lar to the 1995-1997 
results found bv IShara et all (120121) . 

The vario us properties of the re mnant were discussed 
in detail bv lSeaauist et al . (11989a ). The long axis of the 
bipolar H I and IRAS 100 /im map lies in the northwest- 
southeast direction, and essentially orthogonal to the 
northeast-southw est axis (o f symmetry o f the X-ray rem¬ 
nant, see below). iSeaquist et al.l (Il989al ) concluded that 
the nova ejecta are interacting with a cone of this pre¬ 
existing ambient medium in the southwest where the 
remnant appears brightest. This region also corresponds 
to a particularly bright structure lying at an angular sep¬ 
aration of approximately 0'/5 -l' from the ce ntra l star in 
the 19 0 1 and 1902 images bvlRitchevI (|1902j ) and lPerririel 
<|1902f) . ISeaquist et all () 1989 ah suggested this might be 
a cone of material eje cted during the planetary nebu¬ 
lar phase suggested by IBode et all (Il987t) and that the 
remna nt is being dece l erated by interaction with this ma¬ 
terial. ISeaquist et al.l (Il989af) also showed that the rem¬ 
nant expansion was consistent with an adiabatic evolu¬ 
tion in a medium with density declining with the in verse 
square of the radial distance. Here the ISeaquist et al.1 
d!989af) analysi s was based on data obt ained up until the 
late 1980’s, and lAnuoama fe Kantharial (120051) discussed 
subsequent observations made until 2003, presenting ev¬ 
idence that the deceleration of the blast expansion had 
decreased since 1987. They concluded that the density 
into which the southwest quadrant was expanding in the 
later data was lower that in 1987 by a factor of about 7. 

Regarding non-optical wavelengths, non-thermal radio 
emission indicating the presence of accelerated electrons 
interacting with the remnant mag netic field was detected 
in the southwest quadrant by iRevnolds fc Chevalierl 
(1984a), who estimated that a minimum of 1% of the 
10 45 erg explosion energy has gone into the non-thermal 
electrons and the magnetic field. In the X-ray regime, 
GK Per has the only nova remnant for which clumps were 


TABLE 1 

Observation Logs and Spectral Properties ' 1 




Obs.ID 

Par. 

Unit 

650 

15741 

t b 

(d) 

36148.0 

41182.2 

t b 

(yr) 

99.0 

112.8 

fstr C 

(UT) 

2000-02-10 

2013-11-22 



01:15:21 

20:51:14 

^end 

(UT) 

2000-02-11 

04:46:38 

2013-11-24 

00:44:32 

£exp° 

(ks) 

95.3 

97.4 

8o ff d 

(arcmin) 

1.25 

0.18 

SC e -‘ 

(counts) 

3693 

793 

NR e,f 

(counts s -1 ) 

3.06xl0 -2 

5.68xl0 -3 

ME e,f 

HR E 

(keV) 

0.61 

-0.68i0.10 

0.66 

—0.73i0.10 

N h 

k B T 

Zo 

T 

F EM h 

(10 21 cm" 2 ) 

(keV) 

(solar) 

(solar) 

(10 9 s cm -3 ) 

(10 54 cm" 3 ) 

2.10 (fixed) 

r\ 9q+0-03 
U.ZO—0.04 
f) oi+0.09 
u - 01 —0.06 
2.84lj;“ 
4-82 +HI 

o+0.48 

i,1 ®-0.22 

2.10 (fixed) 

pi OO+0.05 

u>zz -0.05 

0 in - * -0 09 

U. l D _o.05 
i n-\ +1.32 
‘ 0.93 

Q 89 

* 5 *' ' -1.57 

i-63±S;gS 

Fx‘ 

L x f ’ h 

(10 -14 ergs s 1 cm -2 ) 

(10 3 ° ergs s 1 ) 

8-09tJ;$ 

9.70tj;?° 

4.847j; 99 

5.76ij;®® 

X^/d.o.f. 

38.4/28 

10.7/14 


“Statistical uncertainties indicate 1 <t confidence ranges. 

“Elapsed days and years from the discovery of the nova (15436.6 d 
in modified Julian date) to the middle of each observation. 

“Start dates (fstr), end dates (t en d); an d exposure times (t e xp). 
d Off-axis angles (0 o u) of the white dwarf binary, of which the 
nova explosion in 1901 caused the surrounding remnant. 

“Source counts (SC), net rates (NR), and median energies (ME). 
^Values are estimated in the 0.5—1.2 keV energy band. 
s Hardness ratios (HR) of the remnant spectra. The values were 
defined as (H—S)/(H+S), where H and S are photon fluxes in the 
harder (0.8-1.2 keV) and softer (0.5—0.8 keV) bands, respectively. 
h Values are for a distance of 470 pc HMcLaughlinl[19601 ) ■ 

spatially resolved. ROSAT discovered the extended neb¬ 
ula for the first time in 1996 (iBalman fc Qgelmanll 19991) . 
finding the bulk of the X-ray emission came from the 
southwestern quadrant, then Chandra revealed the pres¬ 
ence of local X-ray stru ctures in a “ lumpy and asym¬ 
metric nebula” in 2000 dBalmanl 120051) . The Chandra 
ACIS spectrum was dominated by a softer optically- 
thin thermal plasma component below 1.0-1.2 keV, but 
a h arder non-thermal component was also found by 
iBalmanl (120051 ), who concluded that the remnant was still 
in the adiabatic phase in which the radiated energy loss 
is still negligible in comparison with the initial explosion 
energy. 

3. OBSERVATIONS AND REDUCTION 

Chandra has now visited GK Per twice for imag¬ 
ing spectroscopy using the ACIS-S3 back-illuminated 
CCD. The first exposure was taken on 2000 Febru¬ 
ary 10 (Obs.ID = 650) with a total observing time of 
95.3 ks. Our additional observation was obtained on 
2013 November 22 (Obs.ID = 15741) with an exposure 
of 97.4 ks (see Table Q] for details). Both data sets were 
re-processed by the standard pipeline using the Chan- 
dra Interactive Analy sis of Observations (CIAO) package 
(|Fruscione et al.ll2006D with calibration database version 
4.6.2. 
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Fig. 1.— Chandra ACIS-S3 0.5-1.2 keV photon flux images of 
GK Per in (a) 2000 and (b) 2013. Data were smoothed by a Gaus¬ 
sian blur with a 2cr radius of 5 pixels (2'/46), and were rendered 
with a logarithmic intensity greyscale covering up to 2x10“ 7 pho¬ 
tons cm -2 s -1 . The source and background extraction regions are 
indicated by the solid and dashed annuli, respectively. Areas of 
diffuse emission identified by the vtpdetect algorithm are shown 
as green polygons. The 45° blue sectors show the source and back¬ 
ground extraction regions that were rotated in a clockwise direction 
to investigate the azimuthal dependence of the X-ray flux (see text 
and Figure |5|). The red rectangles of size 50 // x20 // parallel to a po¬ 
sition angle of 122° are the regions used for the projection profiles 
of Figure |3] On the lower image, the edge of the CCD chip is also 
shown by the solid black line. 


4. ANALYSIS 

4.1. Diffuse Emission and Comparison 

As found by iBalmanl (120051 ) from the 2000 February 
data, the latest observation reveals the brightest X-ray 
emission to lie in the southwest quadrant and to be char¬ 
acterized by a prominent arc with extended “wings” of 
emission running from the northwest to the southeast. 
The source is notably fainter than it appeared in 2000. 

We limit detailed analysis to soft X-ray events in 
the 0.5-1.2 keV energy band. At lower energies, the 
detector efficiency, especially around the carbon K- 
edge at 0.28 keV, is significantly attenuated by a layer 
of molecular contamination that appears to be de- 
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Fig. 2.— Comparison between radio, optical, and X-ray images 
of GK Per rendered in red, green, and blue, respectively, on loga¬ 
rithmic scales tuned to clarify features of interest. The radio data 
was taken by the B-configuration VLA in 1997 using the L-band 
receiver at 1.45 GHz. The optical data was by the WFPC2 instru¬ 
ment of the HST in 1995 using the Ha+[N II] filter. The X-ray 
image is by Chandra in 2000; the same data as the panel (a) in 
Figure 2] 


positing on the instrument filter (IMarshall et al.l 12004 
\Chandra X-rav Center! I2010D . At higher energies, the 
source fluxes become comparable to, or lower than, the 
background levels. The background-subtracted X-ray 
count rates were estimated to be 3.06±0.07x 10~ 2 and 
5.68±0.32xl0” 3 counts s -1 in 2000 and 2013, respec¬ 
tively, corresponding to a decline of 81±1%. Based on 
folding a model fit to the observed spectrum (see below) 
through the instrument effective areas in the two epochs, 
a drop of 70% in count rate is expected due to the decline 
in the detector efficiency at low energies. We attribute 
the difference between the observed and expected decline 
to the remnant development over the intervening 13.8 
years, indicating an X-ray flux drop of 38±5%. No large 
differences in median photon energy or hardness ratio 
were found between 2000 and 2013 (see also Table □ for 
details). 

Figure Q] shows 0.5- 1.2 keV photon-flux images ob¬ 
tained by Chandra ACIS-S3 in both 2000 and 2013 
epochs. Exposure maps used for flat-held correction 
were computed at 0.6 keV, corresponding to the me¬ 
dian detected photon energy, using the fluximage soft¬ 
ware in CIAO. A slid ing cell searc h for source candidates 
fcelldetect: [Harnden et al.lll984l) was performed, and 
astrometry was fixed by matching the observed X-ray po¬ 
sitions of three p oint sources with Two-Micron All Sky 
Survey (2MASS; iCutri et all 12003 1 catalogue positions 
of their infrared counterparts: J03311201+4354154 (i.e., 
the white dwarf binary of GK Per), J03312219+4356461, 
and J03310820+4357503. The celldetect position es¬ 
timates have precisions better than 0"2, and the result¬ 
ing relative astrometric precision between the two epochs 
based on the three reference sources remains at a similar 
level. The distribution of diffuse sources was determined 
from a Voro noi Tessellatio n and Percolation algorithm 
fvtpdetect: lEbeling fe Wiedenmannlll993l) applied to a 
combined exposure-corrected image in which the two ob- 
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Fig. 3.— Radial projection profiles in the 0.5—1.2 keV band de¬ 
rived from the rectangular regions in Figure Q] smoothed by Gaus¬ 
sian blurs of 5 pixel radius at 2a, are shown color-coded (black for 
2000 and red for 2013). The vertical dotted lines mark the posi¬ 
tions of the peak of the rim emission and the dashed-and-dotted 
lines are the half-maximum levels of each profile. Both indicate 
that the observed expansion of the X-ray nebula is approximately 
1"-2" in 13.8 years. The typical la uncertainties of the data points 
are represented by the vertical bars. 

servations in 2000 and 2013 were merged. 

We then compared the Chandra X-ray data with the 
HST optical image and the Very Large Array ( VLA) 
radio observations in Figure [2] The HST data were ob¬ 
tained from the Hubble Legacy Archive, and we retrieved 
the Ha+[N II] image taken by the Wide Fiel d and Plan¬ 
etary Camera 2 (WFPC2) in 1995 — see IShara et al.l 
( 2012 1 for further details. T he VLA obs e rvat ions 
were discussed by e.g. , Reynolds fo Chevalied (11984H) : 
lAnupama fe Kantharial ( 2005f k Balmanl (12005H . and the 
1.45 GHz L-band image taken by the B-configuration 
in 1997 was obtained from the National Radio Astron¬ 
omy Observato ry VLA Archive Survey (NVAS) system 
(jCrosslev et al.1 120081) . We combined the data from the 
different wavelengt h regi mes using the SAOImage DS9 
image analysis tool ( Jovell20lTh . rendering radio, optical, 
and X-rays in red, green, and blue, respectively. The 
conspicuous coincidence of the radio and X-ray emis¬ 
sion, concentr ated at t he b right rim to the southwest, 
was noted bv iBalman! ( 200 5b who also presented com¬ 
parisons with optical image contours. 

4.2. Limits on Hard X-ray Evolution 

Unfortunately, our additional data provide no further 
insights i nto t he n on-the rmal hard X-ray signatures re¬ 
ported bv IBalman! (2 0051 ) because of the background con¬ 
tamination above about 1 keV. We investigated this non- 
thermal emission by first fitting the extracted spectra 
from both epochs in the 1.2-7.0 keV range using an ab¬ 
sorbed power-law model and estimating the la statisti¬ 
cal uncertainties in the power-law normalization. Using 
these, we calculated count rates in the 0.5-1.2 keV range 
by extrapolation. The fitted power-law slopes were ap¬ 
proximately zero, and the resulting upper limits to the 
contribution of power-law component to the count rate 
are 10% and 23% for the first and second observations, 
respectively. The larger limit from the second epoch re¬ 
sults from the larger uncertainty in the fitted slope for 
those data. However, it is unlikely that the harder com¬ 
ponent has increased in intensity or hardness while the 
rest of the remnant has faded significantly. To place a 
best estimate on the soft X-ray power-law emission con¬ 
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Fig. 4.— Comparison of background-subtracted SB radial profiles 
of GK Per in 2000 (black) and 2013 (red) in the 0.5—1.2 keV band 
computed from the exposure-corrected photon-flux images, taking 
into account the degradation of the quantum efficiency of ACIS-S3. 
The dashed lines show 0.6 keV PSFs that were normalized to the 
maxima of the observed peaks, and further to the minima of the 
background levels. The lower panel indicates the ratios of 2013 to 
2000 data, suggesting that the SB drop was 36±6% in the 20”- 
40” plateau area, which is indicated by the orange region and the 
dashed-and-dotted lines. 

tribution, we therefore imposed additional constraints of 
not allowing either the spectrum to get harder or the nor¬ 
malization to increase from the first epoch to the second. 
We introduced a lower limit to the power-law slope and 
an upper limit to the normalization corresponding to the 
ler limits derived from the first observation and found a 
lcr limit of an 18% contribution due to the harder com¬ 
ponent, although the actual contribution is likely similar 
to the 10% one of the earlier epoch. Since the power-law 
spectrum would be essentially flat, it has little influence 
on the interpretation of the soft X-ray data and we dis¬ 
regard it in the remainder of the analysis. 

4.3. Radial Profiles and Remnant Expansion 

In order to examine the expansion of the remnant be¬ 
tween the two Chandra observation epochs, the projected 
and summed profile of the bright southwest rim was an¬ 
alyzed. A rectangular extraction region was placed on 
top of, and tangential to, the rim (indicated by the red 
rectangles lying to the southwest in Figure [J), and X-ray 
events were extracted and summed in the tangential di¬ 
rection such that the resulting histogram shows the vari¬ 
ation of integrated flux (IF) in an approximately radial 
direction from the nebula center. 

Figure [3] illustrates the projection profiles of these 
brightest regions in the X-ray images of both epochs. 
Comparing the profiles at their peaks suggests an ob¬ 
served expansion of 1"5, while at half maxima linear in¬ 
terpolation indicates 1"9 in 13.8 years. These two esti¬ 
mates are significantly larger than any systematic error in 
the relative positional accuracy of the two observations. 
In the following, we adopt the measurement based on the 
half maxima for discussing the remnant expansion. 

Radial surface brightness (SB) profiles derived by az- 
imuthally summing the observed signal as a function of 
radial distance from the central binary position, with a 
0'/5 step size, are shown in Figure [4] Also shown is a 
0.6 keV point spread funct ion (PSF) calc ulated using the 
Chandra Ray Tracer tool dCarter et al.ll2003D . Compar¬ 
ison of the central source counts as a function of radius 
with the PSF prediction indicates the source is heavily 
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Fig. 5.— Comparison of 0.5-1.2 keV background-subtracted az¬ 
imuthal dependence in integrated flux in 2000 (black) and 2013 
(red). The solid and dash-dotted lines are the results of using dif¬ 
ferent background estimates; the former used the one-eighth out¬ 
ward sectors at each azimuth angle, of which the background levels 
are shown by the dashed lines, while the latter is from the entire 
background annulus (see Figure [l] for details). The lcr confidence 
ranges of the former estimates are illustrated by the color-coded 
bands. The lower panel indicates the ratios of 2013 to 2000 data. 

piled up in the center. The excess counts over the PSF 
model out to about 10" are then expected to arise as a 
result of the PSF being normalized to the peak counts, 
although it is possible some of the excess could be ex¬ 
tended emission from the remnant. 

Plateaus in the SB 20"-40" from the center indicate 
regions where the nebula emission clearly dominates any 
signal from the central source. Fitting a constant SB 
level to these segments of the radial profiles resulted in 
best fits of 1.48±0.04x 10 -8 and 0.95±0.06xl0 -8 pho¬ 
tons cm -2 s -1 arcsec -2 in 2000 and 2013, respectively, 
corresponding to a decay of 36±6% in 13.8 years. The ra¬ 
dial profiles decay to the background levels at distances 
of 50"-60". Based on the radial profiles, the remnant 
signals for further spectral analysis were extracted from 
annular regions between 20"-60", avoiding potential con¬ 
tamination from the PSF wings. Background events were 
accumulated from the 70"-100" annular region, exclud¬ 
ing the areas where diffuse emission was detected by the 
vtpdetect algorithm (see Figure |T] for details), and the 
background rate was computed by normalizing to the ra¬ 
tio of net source and background extraction regions. 

4.4. Azimuthal Dependence and Evolution 

The azimuthal dependence of the nebular emission 
was investigated by assessing the source and background 
counts in 45° segments of their respective annular extrac¬ 
tion regions and rotating these segments in the clockwise 
direction with a 1° step size. The resulting IFs as a 
function of azimuth computed using both the rotating 
background sector and just the average background in 
the annulus (again with diffuse emission excluded) are 
illustrated in Figure [5j The results reflect the “by-eye” 
assessment that the brightest part of the remnant is in 
the south and the faintest in the east. The ratios between 
2000 and 2013 curves indicate that the flux decline is in¬ 
versely related to X-ray brightness - fainter regions are 
getting more faint at a faster rate than brighter regions. 
A simple average of the IF ratios is about 70%, with a 
maximum ratio of 90% in the bright area to the south, 
and down to 20% to the northwest. 


Fig. 6.— Background-subtracted spectra of the GK Per remnant 
in 2000 (black) and 2013 (red). The best-fit models of each spec¬ 
trum are shown by the dashed lines. The lower panel shows the 
residuals from the best-fit models. The energies of plausible Kq 
emission lines are labeled with solid and solid-and-dashed lines for 
H- and He-like ions, respectively. 

4.5. Spectral Analysis 

As no significant secular trends, except for background 
variations, were found in the light curves constructed 
from the background-subtracted remnant events, spec¬ 
tral analysis proceeded based on the whole of each data 
set. Background-subtracted X-ray spectra from both 
epochs are illustrated in Figure [G] Emission from Ne IX 
Ka was clearly present in both epochs, while O VII 
emission is also possibly discernible. We applied a fit¬ 
ting model comprising emission from a non-equilibrium 
ionization collision-dominated plasma plus photoelectric 
int erst ellar a bsorp tion (Tuebingen-Boulder absorption; 
iWilms et ali l2000i ) with a hydrogen-equivalent column 
density (Ah) fixed at 2.1xl0 21 cm -2 . The absorption 
column was not well-constrained by the Chandra data 
and the latter was estimated from the H I maps of 
iKalberla et al.1 (120051) for a 1° cone radius around the 
GK Per position. Here we experimented first with al¬ 
lowing the plasma abundances of C, N, O, and Ne to 
vary from their solar values, as might be expected in the 
shell of a planetary nebula and also in the ejecta of a 
classical nova. As a result no significant constraint on 
the C abundance has been given from the spectral fit¬ 
ting. Regarding the N emission, the best-fits and lcr 
statistical errors suggest an abundance of 0.3-2.7 times 
the solar value. The reason for such poor constraints 
is that the energies of the Ka emission lines from those 
metals are slightly below the useful low energy cut-off of 
the X-ray spectra. We confirmed that varying the C and 
N abundances provides no significant difference on other 
parameters, and thus the elemental abundances of the 
emission component were fixed to be solar, except for O 
and Ne from which the energies of the Ka lines are in 
the range of the extracted spectra. For the final spec¬ 
tral fitting, the free parameters were then: (1) plasma 
temperature (ksT), (2) elemental abundances of O (Zq) 
and Ne (Zx e ), (3) the non-equilibrium plasma ionization 
timescale (r), and (4) a volume emission measure (Vem)- 
The observed fluxes (Fx) and absorption-corrected lu¬ 
minosities (Lx) in the 0.5-1.2 keV energy band were 
also determined. The flux decay between the 2000 and 
2013 epochs determined from these spectra amounts to 
40%. The spectral properties and the best-fit parame¬ 
ters producing minima in the y 2 test statistic defined by 
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iGehrelsI (|1986T ) are summarized in Tabic [1] The results 
for the 2000 epo ch are broadly consistent with the re¬ 
sults of IBalmanl (120051) . and the fitted temperature for 
both epochs of 0.2 keV is similar to the 2xl0 6 K tem¬ 
perature derived from the 1992 ROSAT observation by 
iBalman fe Ogelmani ( 1999 ). Independent spectral anal¬ 
ysis of different regions of the nebular was investigated 
but not pursued due to insufficient signal-to-noise ratio. 


5. DISCUSSION 

5.1. Energetics 

We start with an overview of the energetics of the ex¬ 
plosion that lead to the GK Per remnant. First, assuming 
the angular radius R = 1' corresponds to 0.14 pc at a dis¬ 
tance of 470 pc, the upper limit to the remnant volume 
assuming a spherical shape is V = 3.4xl0 53 cm 3 . The 
geometry of the emission is instead presumably a thin 
shell with thickness R/ 12, as expected for a compression 
ratio of 4 for a strong shock, which suggests a plasma 
emitting volume of V = 7.8xlO 52 cm 3 . Here, the ob¬ 
served volume emission measures, Hem = n e n\V, where 
n e and are electron and ion number densities, respec¬ 
tively, are 1-2xlO 54 cm -3 at the distance 470 pc. The 
larger uncertainties on the fitted Vem values result from 
degeneracy with the weakly-constrained Ne and O abun¬ 
dances. Further assuming solar abundances, the electron 
to ion ratio is n e = 1.18ni, and so rij ~ 4 cm” 3 , which 
corresponds to the pre-shock total number density of no 
~ rii/4 ~ 1 cm -3 . The swept-up mass is M ~ m^n\V 
~ 2xl0~ 4 Mq, where toh is a hydrogen mass. _ 

If the remnant is in an adiabatic stage ( Scaquis t et all 

Il989al IBalmanl 120051 : lAnupama fe Kantharial I2005D . in 

which radiating energy loss is still negligible compared 
with the initial explosion, the plasma state can be com- 
pared with the Sedov-Taylor solution (iSedovI fl959l ) rep¬ 
resented by the following three equations, 




10 4 yr v 10 51 ergs' v lcm -3 ' 


s cm 


_l. cm. 


„ = 5x10 7 (^-)”5(-^- ) 5 ( n ° I »I 

v 10 4 yr 3 3 10 51 ergs 3 v lcm -3 s 

T = 3 x 10 6 ( —)“®( E 


(1) 

( 2 ) 


( 3 ) 


10 4 yr * * * * v 10 51 ergs' v lcm _3 


Here, the simultaneous equations consist of the six 
parameters: the radial distance to the outer shock 
R, the blast-wave velocity v, the shock temperature 
T, the elapsed time t, the total explosion energy E, 
and the total number density no. The velocity is v 
= (2/5 )(R/t) ~ 500 km s -4 with R = V and t = 
112.8 years, and similarly the temperature is k^T = 
1.8xl0 5 (i?[pc]/t[yr]) 2 ~ 0.28 keV. The radial distance is 

then R = 0.31(.E[10 51 ergs]/n 0 ) 1//5 t[yr] 2 ^ 5 pc, which can 
be converted to the explosion energy E ~ 2xl0 45 ergs 
using no = 1 cm” 3 . This distance is fairly insensi¬ 
tive to the exact particle density assumed, The ob¬ 
served electron temperature of fee T ~ 0.2-0.3 keV is 
almost consistent with the Sedov-Taylor estimates, sug¬ 
gesting that the electrons and ions are near equilibration. 
All these estimates lie within typical ranges for nova 
explosion s obtained from both observations and the¬ 
ory (e.g.. iPrialnik fe Kovetd 1 199.51 iStarrfield et all 120091 : 
iMetzeer et al.ll2014D . 


5.2. Proper Motion 

The observed expansion of the X-ray brightest part of 
the remnant in the 13.8 years separating the Chandra 
observations derived from the radial-profile comparisons 
is l'/9, corresponding to an expansion rate of 0'/14 yr _1 . 
Assuming a distance of 470 pc, the velocity of the out¬ 
ermost X-ray emitting region is 300 km s _1 . This is 
slightly lower than that the 600-1000 km_s -1 of o ptical 
knots moving in a similar direction (IShara et all 120121 : 
iLiimets et al .1120 lH . see also lDuerbecklll987D . as well as 
the Sedov-Taylor estimate of 500 km s _1 . The veloc¬ 
ity inconsistencies might imply that the shock speed was 
faster than the estimate from the apparent X-ray advance 
and that this measurement is not representative of the 
advance of the shock front itself. Further, as we discuss 
in Sect. 15.31 the expansion of the blast wave must have 
been non-uniform. 

5.3. Cooling or Thinning? 

The flux estimates based on radial profiles, the IF az¬ 
imuthal dependence, and the spectral fitting all indicate 
that the 0.5-1.2 keV brightness has been reduced on av¬ 
erage to 60-70% of its value in 13.8 years, corresponding 
to a decline rate of about 3% yr _1 . This is very similar to 
the flux decay rate of the opti cal knots in the remnant of 
2.6% yr -1 (ILiimets et al.ll2012l) . and to the annual secu¬ 
lar decrease of 2.1% in the flux density at 1.4 and 4.9 GHz 
betwe en 1984 and 1997 found bv lAnunama fe Kantharial 
(I2005D . There is no obvious significant X-ray spectral 
evolution accompanying this decay, and all the fitted 
spectral model parameters for the two epochs were sta¬ 
tistically consistent (see Table [I]). 

For an emitting plasma with cosmic abundances, the 
observed temperature of 0.2-0.3 keV implies a total 
coolin g coefficient of orde r 10~ 22 -10” 23 erg cm 3 s -1 
(e.g., iGehrels fe Williams! I1993D . Here, the dorni- 
n ant cool in g occurs from the line emission of Fe 
dBoehringer fc Henslerl 119891 ). Since the shocked cir- 
cumstellar environment is likely to have a chemical 
composition more like that of typical planetary neb¬ 
ulae, the actual cooling coefficient could differ from 
this value — the tendency of the Fe depleti on in 
galactic pla neta ry nebula e (iDelgado Inglada et al.ll2009fc 
iDelgado-Inglada fc RodriguezI 120141) suggests that the 
plasma cooling becomes less effective also in the GK Per 
remnant. The cooling timescale is therefore orders of 
magnitude longer than the present lifetime of the rem¬ 
nant. As the temperature was almost constant dur¬ 
ing 2000-2013 as expected, the radiated energy inte¬ 
grated over 13.8 years is of order 10 39 10 41 ergs or lower. 
This strongly suggests that radiative cooling is negli¬ 
gible in comparison with the explosion energy of order 
10 45 ergs, and supports earlier conclusio ns t hat th e rem ¬ 
nant is in the adiabatic phase (ISea.auist et al.l [l989a; 
lAnupama fe Kanthariall2005l : lBalmanll2005l) . 

We thus conclude that the most plausible cause of the 
30-40% decline in emission between 2000-2013 is not 
the plasma cooling but just the thinning and expan¬ 
sion of the emitting nebula. Assuming the blast veloc¬ 
ity of 500 km s -1 , the emitting volume increased 20% in 
13.8 years. Were the blast wave expanding into a medium 
of uniform density, the increase in swept-up mass com¬ 
bined with an increase in the shocked gas volume in the 
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expression for the emission measure, means the X-ray 
emission should increase with the cube of the remnant 
radius. In the Sedov phase, since R oc t 2,/5 , the emission 
measure should increase with time as Vem oc f 6 / 5 . That 
this is definitely not the case indicates that the medium 
into which the GK Per remnant has been recently ex¬ 
panding must be non-uniform. Indeed, as discussed in 
Sect.[2l the complex non-uniform nature of the circumbi- 
nary medium of GK Per has been evident for more than 
a century. The observed decline in X-ray emission cor¬ 
responds best to the case of negligible swept-up mass 
during 2000-2013 and therefore a number density at the 
shock front in the later epoch much lower than at earlier 
times. The later volume emission measure Vem = n e n\V 
would then be about 80% of the 2000 value. 

The inverse relation between the X-ray brightness and 
the decline rate found in the azimuthally-dependent data 
indicates that the currently fainter regions are fainter 
because they have dimmed more rapidly, and continue 
to do so. This is especially prominent to the east, at an 
azimuthal angle of around 270°. These regions might be 
fading faster because they are expanding more rapidly, 
although the X-ray data are insufficient to test this. 

Asymmet ric expa nsion is not see in the optical knots 
(lLiimets et al.lf2012 ) which actually experienced a recent 
brightening in the east. Such divergent behaviors are not 
unrealistic since optical knots are probably dominated by 
shocked ejecta and we expect the X-ray emission to in¬ 
stead originate largely from the shocked medium. Com¬ 
parison of the X-ray and optical emission in the multi¬ 
wavelength image in Figure [2] reveals that some optical 
knots precede the X-ray emitting shock front. Similar 
behavior has been seen in optical knots of the Cas A 
supernova remnant (lHammell fe Fesenll2008f) . The main 
shell of the brightest knots is quasi-spherical, and close 
inspection of the X-ray emission indicates that it is co- 
spatial with knots. This indicates that the knots are 
likely not the same material as the X-ray emitting plasma 
that has condensed and cooled, and is consistent with 
the knots being shocked ejecta (e.g. , iBode et all 12004 
lLiimets et al.ll2012t iShara et all 1 20121) . 

The most obvious departure of optical and X-ray mor¬ 
phology is in the X-ray “wings” pointing to the north and 
east that extend well beyond the optical knots. IBatmanl 
(I2005T) noted that the expansion rate for these wing re¬ 
gions must be about 2600-2800 km s -1 - more than 
twice the rate of the optical knots and the bright emis¬ 
sion in the southwest — and drew comparison to the 
“cone” of earlier ejected material propos e d in th e plan¬ 
etary nebula scenario by iSeaauist et all (l!989aT) . The 
implication would be that we are still seeing largely the 
shocked medium, but that the density experienced by the 
blast in these wing directions up until encountering the 
denser material must have been much lower than in the 
southwest direction to permit the more rapid expansion. 

6. CONCLUSIONS 

Comparison of Chandra X-ray imaging spectroscopy 
observations of the remnant of the GK Per classical nova 
event of 1901 obtained in 2000 and 2013 has revealed the 
following results. 

1. The remnant faded significantly in the interval be¬ 
tween the two observations, with a mean flux de¬ 


cline of 30-40%. The rate of fading is spatially 
dependent, with more rapid fading generally coin¬ 
ciding with regions of lower surface brightness. 

2. Comparison of X-ray and optical data confirms 
qualitatively the very different morphology of cool 
knots and hot gas. The highly asymmetric X-ray 
remnant contrasts with a quasi-spherical shell of 
complex knots. Some optical knots appear to be 
slightly ahead of the X-ray emitting shock front, 
and also appear X-ray dark. 

3. The expansion of the remnant deduced from the 
difference in the bright rim of X-ray mission to the 
southwest in the two different epochs is l'/9, or an 
expansion rate of 0"14 yr -1 . This corresponds to 
velocity of approximately 300 km s -1 for a distance 
to GK Per of 470 pc. 

4. Assuming the X-ray emission arises from a rela¬ 
tively thin shell of shocked plasma, the swept-up 
mass is approximately 2xl0~ 4 M 0 and the ion 
density is 4 cm -3 . A Sedov-Taylor solution yields 
an explosion energy of 2xl0 45 ergs, and electrons 
and ions are close to temperature equilibration. 

5. Since the radiative cooling time of the shock-heated 
gas is much longer than the age of the remnant, 
we conclude that the X-ray fading is a result of 
adiabatic expansion at a rate of about 3% per year 
in the southwest rim. The more drastic fading in 
regions of low surface brightness, such as in the east 
of the remnant, is then expected to be caused by 
more rapid expansion into a lower-density medium. 

6. The general evolution of the remnant is qualita¬ 
tively consistent with expectations based on early 
post-explosion photographs of pre-existing nebu¬ 
losity, with a highly-structured and inhomogeneous 
shell-like morphology. The fading observed be¬ 
tween 2000 and 2013 is then caused by the re¬ 
cent expansion of the remnant into regions of much 
lower density. This is consistent with a planetary 
nebula origin for the pre-blast circumstellar mate¬ 
rial illuminated by light echoes from the explosion 
itself. 


We thank Chandra for allocating a part of the tele¬ 
scope time in Cycle-15, and acknowledge HST and 
VLA for their archival images. D.T. thanks T.Yuasa, 
K.Morihana, and J.Ueda for valuable discussions, and 
acknowledges financial support from the RIKEN/SPDR 
program, Grant-in-Aid for JSPS Fellows for Research 
Abroad, and Chandra grant GO4-15025X. J.J.D. and 
P. S. were supported by NASA contract NAS8-03060 
to the Chandra X-ray Center, and thank the Director, 
B. Wilkes, for continuing advice and support. 



















D. Takei et al. 


REFERENCES 


Anupama, G. C. & Kantharia, N. G. 2005, A&A, 435, 167 
Anupama, G. C. & Prabhu, T. P. 1993, MNRAS, 263, 335 
Balman, §. 2005, ApJ, 627, 933 
Balman, §. Sz Ogelman, H. B. 1999, ApJ, 518, Llll 
Bode, M. F., O’Brien, T. J., & Simpson, M. 2004, ApJ, 600, L63 
Bode, M. F., Roberts, J. A., Whittet, D. C. B., Seaquist, E. R., & 
Frail, D. A. 1987, Nature, 329, 519 
Boehringer, H. & Hensler, G. 1989, A&A, 215, 147 
Carter, C., Karovska, M., Jerius, D., Glotfelty, K., & Beikman, S. 
2003, in Astronomical Society of the Pacific Conference Series, 
Vol. 295, Astronomical Data Analysis Software and Systems 
XII, ed. H. E. Payne, R. I. Jedrzejewski, & R. N. Hook, 477 
Chandra X-ray Center. 2010, Chandra Technical Information 
Crossley, J. H., Sjouwerman, L. O., Fomalont, E. B., Sz Radziwill, 
N. M. 2008, in Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, Vol. 7016, Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, 0 

Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR 
Online Data Catalog, 2246, 0 
Delgado-Inglada, G. & Rodriguez, M. 2014, ApJ, 784, 173 
Delgado Inglada, G., Rodriguez, M., Mampaso, A., Sz Viironen, 

K. 2009, ApJ, 694, 1335 
Duerbeck, H. W. 1987, Ap&SS, 131, 461 
Ebeling, H. Sz Wiedenmann, G. 1993, Phys. Rev. E, 47, 704 
Fesen, R. A., Hammell, M. C., Morse, J., et al. 2006, ApJ, 645, 
283 

Fruscione, A., McDowell, J. C., Allen, G. E., et al. 2006, in SPIE 
Conference Series, Vol. 6270, SPIE Conference Series 
Gehrels, N. 1986, ApJ, 303, 336 
Gehrels, N. Sz Williams, E. D. 1993, ApJ, 418, L25 
Hammell, M. C. Sz Fesen, R. A. 2008, ApJS, 179, 195 
Harnden, Jr., F. R., Fabricant, D. G., Harris, D. E., Sz Schwarz, 

J. 1984, SAO Special Report, 393 
Joye, W. 2011, in Astronomical Society of the Pacific Conference 
Series, Vol. 442, Astronomical Data Analysis Software and 
Systems XX, ed. I. N. Evans, A. Accomazzi, D. J. Mink, Sz 
A. H. Rots, 633 

Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, 
A&A, 440, 775 


Katsuda, S., Long, K. S., Petre, R., et al. 2013, ApJ, 763, 85 
Lawrence, S. S., MacAlpine, G. M., Uomoto, A., et al. 1995, AJ, 
109, 2635 

Liimets, T., Corradi, R. L. M., Santander-Garcia, M., et al. 2012, 
ApJ, 761, 34 

Marshall, H. L., Tennant, A., Grant, C. E., et al. 2004, in Society 
of Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, Vol. 5165, X-Ray and Gamma-Ray Instrumentation for 
Astronomy XIII, ed. K. A. Flanagan & O. H. W. Siegmund, 
497-508 

McLaughlin, D. B. 1960, in Stellar Atmospheres, ed. 

J. L. Greenstein, 585 

Metzger, B. D., Hascoet, R., Vurm, I., et al. 2014, MNRAS, 442, 
713 

Patnaude, D. J. & Fesen, R. A. 2009, ApJ, 697, 535 

Perrine, C. D. 1902, PASP, 14, 151 

Prialnik, D. Sz Kovetz, A. 1995, ApJ, 445, 789 

Reynolds, S. P. Sz Chevalier, R. A. 1984a, ApJ, 281, L33 

Reynolds, S. P. Sz Chevalier, R. A. 1984b, ApJ, 281, L33 

Ritchey, G. W. 1902, ApJ, 15, 129 

Schaefer, B. E. 1988, ApJ, 327, 347 

Seaquist, E. R., Bode, M. F., Frail, D. A., et al. 1989a, ApJ, 344, 
805 

Seaquist, E. R., Bode, M. F., Frail, D. A., et al. 1989b, ApJ, 344, 
805 

Sedov, L. I. 1959, Similarity and Dimensional Methods in 
Mechanics 

Shara, M. M., Zurek, D., De Marco, O., et al. 2012, AJ, 143, 143 
Slane, P., Gaensler, B. M., Dame, T. M., et al. 1999, ApJ, 525, 
357 

Slavin, A. J., O’Brien, T. J., Sz Dunlop, J. S. 1995, MNRAS, 276, 
353 

Starrfield, S., Iliadis, C., Hix, W. R., Timmes, F. X., Sz Sparks, 
W. M. 2009, ApJ, 692, 1532 

Uchiyama, Y., Aharonian, F. A., Tanaka, T., Takahashi, T., & 
Maeda, Y. 2007, Nature, 449, 576 
Williams, A. S. 1901, MNRAS, 61, 337 
Wilms, J., Allen, A., Sz McCray, R. 2000, ApJ, 542, 914 
Yamaguchi, H., Ozawa, M., Koyama, K., et al. 2009, ApJ, 705, L6 



